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Abstract 

To understand the mechanism of the fermion pair and fermion-antifermion pair condensation, the 
solutions of Bethe-Salpeter equation in QED3 is examined. In the ladder appoximation our solution 
for the axial-scalar is consistent with Ward-Takahashi-identity for the axial vector currents. Since 
the massless scalar-vector sector is described by a coupled integral equation, it is difficult to solve 
explicitly. We approximate the equation for large and small momentum region separately and con- 
vert them into differential equations in position space. These equation can be solved easily. Boundary 
condition at the origin leads the eigenvalue for dimensionless coupling constant A = e^/m. There ex- 
ists solutions for massless scalar-vector fermion-antifermion (fa) system with discrete spectrum. In 
our approximation massless-scalar- vector ff systemes does not seem to exist. 



101 



Contents 



I. Introduction 102 

II. Spinor-Spinor BS equation 103 

A. Massless boson in QED 103 

B. Massless boson in QED3 106 

III. Ward-Takahashi-identity for Axialvector currents 110 

IV. Summciry 113 
V. References 113 



I. INTRODUCTION 

In the theory of superconductivity and superfluidity BCS model is very familiar and 
useful to analyse their properties [1]. If we assume the S^^^x — y) function type intearction 
between fermions it is not difficult to solve the gap equation and its dependence of the 
coupling const ant. However we have not been known the reason why the electron pair form 
a bound state. On the other hand there exisists solutions of Bethe-Salpeter equation in the 
ladder approximation for fermion-fermion pair with discrete spectrum [2] .Recently a bound 
state for quasi partices(exiton) are considered in terms of approximate Bethe-Salpeter(BS) 
equation in QED3 for phase flucuating d-wave superconductor near the node to measure the 
resonant spin response[3].Thc Schrodinger type equation for particle-hole boundstate with 
potential 1/r^ was derived and the eigenvalue condition emerged. Its solution indicates the 
exsitence of strong coupling phase but seems to be not normalizable.Therefore it is interesting 
to apply relativistic BS equation for fermion pair in the same approximation in[2].We think 
that it is important to solve the equation for the massless boundstate which signals the 
instability of the vacuum under condensation of these bosons in field thoretical model [4]. In 
this work first we examine the existence of solutions for the massless boundstates which are 
normalizable. Since scalar- vector sector is written by coupled integral equation, it is difficult 
to solve the equation explicitly. Therefore we approximate the equation for the large and 
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small momentum region separately. These equations are easily solved in position space and 
we have a correct short and long-distance behaviour. For fa system the baundary condition 
at the origin leads dimensionless coupling constant -\/2e^/47rm = integer where m a fermion 
mass. In section II we introduce spinor-spinor BS equation and show their solutions in QED 
and QEDspj.In section III we compare the solutions of Dyson- Schwinger equation for the 
fermion propagator with axialscalar solutions of BS equation.In our approximation they 
obey the same equation and Ward-Takahashi-identity for axial currents is satisfied. Section 
IV is devoted for summary. 

II. SPINOR-SPINOR BS EQUATION 

A. Massless boson in QED 

BS amplitude in four-dimension is defined [6, 7, 8] 

x{xiX2:B) = {0\T{ij{x,)i^{x,))\B) , (1) 
for the fermion-antifermion bound state \B) — \fa, P^j) with total four momentum or as 

x{xiX2:B) = (0\T{'iP{xi)if'{x2)\B] 



for the fcrmion-fcrmion bound state \B) = |//, P^), where ip^ stands for the charge conju- 
gated field of V'. Homogeneous Bethe-Salpeter equation for fermions in the ladder approxi- 
mation is written in the following form 

x{xiX2:B) = -(? J d^xad'^x^Spixi - X3)x{x3X4::B)Sf{x2 - Xij'y/^D'p' {x3 - x^)^^. (2) 

We can also write the BS equation in a differential form by applying Dirac operator; 

{idi ■ 7 - m)xixiX2:B){id2 • 7 - m) = -e^D'^^{xi - X2)j^ixix 1X2: B)-fi,, (3) 

where m is a dynamical mass and Dp is a photon propagator in the covariant gauge 



D'/iP)-J-- + i^-^)^- (4) 



In momentum space we transform to center of mass and relative coordinate 



X = (Xi X2)/2, X = Xi - X2, (5) 
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^ = Pi+P2,P = (pi-P2)/2. (6) 
Then we define x(P,p),the Fourier transform of tlie Feynman amplitude, by 

x{xiX2:B) ^ exp{iP ■ X) J d'^pexp{ip • x)x{P,p)- (7) 

The BS equation in momentum space assume the form 

P P /" 

for fa system. Since under charge conjugation vector particle is odd, A should be replaced 
by —A. For scalar case the x^{P^P) is written in general which corresponds to spin singlet 
S{P,p) and triplet V{P,p) and tensor T{P,p) 

X'{P,P) = SiP,p)I + P ■ ^V\P,p) +p- ^V\P,p) + a,,T{P,p){P^p, - P^p,). (9) 

For total momentum = case the BS amphtdes are decoupled to scalar- vector and tensor 
and we have a coupled equation for scalar- vector [2, 6, 7, 8]. If we substitute eq (9) to eq (8) 
we get a scalar- vector equation for Euclidean momentum p 

(m^ - - 2pV^(/) = \' j d'q^^„ (10) 

and 

(m^ - P>,V\p') + 2p,S{p') = A^ I ^'^^r^- (11) 

A^ - (3-OA... for ff-system 

= -(3 - OA., for fa-system, (12) 

A^ = 2^ A. .for ff-system ( 13) 

= — 24A..for fa-system, (14) 

where A = e^/(47r)^.ln the Landau gauge {C, = 0) A^ vanishes. The solution is given in 
ref[2].Then eq (11) reduces to a algebraic equation and V^{p'^) is obtained 

V'ip') = , f ,. Sip'). (15) 

From eq (10) and (15) we obtain 

(•"'+P^)-S(p^) = A»/A7^. (16) 



m? — p^ J {p — qy 
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By using the formula(Klein-Gordoii equation for photon) 



□,^-4-^ = -47r^5(^)(p-,) (17) 
we get a differential equation for S{p^) from eq (16) 

where s = p^.The solution is characterized by hypergeometric function 

S{s) = i^^As + m2)-i--/2F(a, A 7; -^), 
[s + m^j m"^ + s 

2 + Vl + 4A^ - Vl + 8A^ ^ 2 + Vl + 4A^ + Vl + 8A^ 
2 '^^ 2 ' 

7 = 1 + Vl + 4A^, (19) 
which satisfy the boundary condition 

H^'^l^SMl^,. (20) 

After angular integration eq (16) is rewritten 
(m^ + s)^ 



m? — s 



S{s) = ^ / dsV>S(s') + ^ / ds'S{s'). (21) 



If we differentiate the above integral equation we get the boundary conditions. In the case 
of a{P) — —n, F is a n-th degree of hypergeometric series in s.We have a discrete set of 
spectrum if A'^ > O(fermion-fermion system) and the eigenvalues are given 

Af = {n + l)[3(n + 1) + VSn^ + 16n + 9], n = 0, 1, 2... (22) 

For the lowest eigenvalue n = (Aq = Ao/3 = 2),wc have 

^0 W ^ — / 2^ 2^5 — ■ 23 
For fermion-antifermion system there exists contineous spectrum for— 1/8<A<0 

= T^^^^' + ^')"'""/'^(«, /?; « + /3 - 7 + 1; ^^)- (24) 
Here we show the profile of Xg{'^) FIG I.Fourier transformation of x{p) is defined 

^^""^ " 4^ / P'^^P^^^^^P^^ 
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for m = 1. In ref[4] introducing auxiliriary field for comosite operator effective action has 
been made. If we minimize it we can derive the Dyson- Schwinger equation for fermion and 
Bethe-Salpeter equation for boundstate. However for the first approximation to the lowest 
excitation, vector was neglected. Thus the BS equation for scalar meson is not a coupled 
equation with vector meson as eq(10),(ll). 

B. Massless boson in QED3 

In ref[3] three dimensional QED is considered as an effective model which discribes the 
two dimensional superconductor, where gap has nodes and the low-energy fermionic excita- 
tions have linear dispersion. The fluctuating phase of superonductor is considered as Berry 
gauge fields which couples to spin degree of freedom of exciton.For these systems we assume 
the validity of relativistic treatment for fermion. For definiteness we use four-dmensional rep- 
resentation of 7 matrix [5]. In ref[3] following type of equation for the eigenvalue e„(|p| ,m) 
and normalized eigenfunction '?/'„(r, |p| , m) for electron-hole system near the nodes were 
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dicsussed 

[-d-d- Veff]ljjn = e„Vn, (26) 
X(P) = ^ I dhe^-ip ■ r) (V^,(r)V^,(r)V^^.(0)V^,(0)> , (27) 

ij 

,)^jy^J«O^M^_ (28) 

^ 4e„(|p| ,m) 

where Vg// is a potential by dressed gauge boson and has an invers-square form for small r 
and x{p) is a sum of boundstate propagator with residue iV'nl^of each pole.The conclusion 
is that there are infinite numbers of negative energy. Here we simpliy notice the results 

i-d-d- = A = (2 + 09^ (29) 

l( ^ ^V43A/4(^/2) 

nz) = 7= . (30) 

It is said that the continuum spectrum that is normalizable into a single boundstate would 

be a so called conformal tower 

rK^Ko{K\z\) 

where —k^ is the renormalized bound-state energy,and Kq is the modified Bessel function 
of the second kind. Apart from these realistic application to condensed matter physics it is 
interesting to solve the massless boundstate problems in three dimension. Now we return to 
relativistic BS equation for fa system to study the same problems. For = 0,T does not 
couple to S, \^.We consider the equation as (16) in the previous section. 

("''+"'''s(p^) = A»/A7^. (32) 



— J {p — q) 

However we cannot solve the above equation as in four dimension since the identity is 

= -47r5(^)(p-?) (33) 

\p-q\ 

in this case.Therefore we expand the left hand side of the eq (32) in p for large and small 
cases. We get 

{-2m' - p')S{p) ^X'J d'qj^^„ (m^ « p'), (34) 
(m^ + 2p')Sip) = X' j d'Q-^^, ip' « m'), (35) 
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where A'^ = e^/27r. We convert these equations in position space 



+ D,)Ss{x) = X''^^, (m |x| < 1), (36) 



\x\ 

m'' -2n^)S l{x) =X^ ^^,{1-^ 7711x1). (37) 



\x\ 

First we solve the following equation for scalar- vector ff system S{x) 

\X\ ZTT 

^ (ir^ r cir 

The relation among S{x), Ss{x), Sl{x) may be clear by scaling the mass or coupling con- 
stant. For simphcity we set m — l.For the ground state I — 0,we have 

^ + -^-S(r)-yS(r) = 0. (39) 

For large r if we neglect terms which are proportional to l/r,we obtain 

(fSjr) 
(Pr 

from this S{r) behaves as e"~^.For small r 



S{r), (40) 



S{r) — u{r)/r. 

then the eq (39) becomes 

|^ + M->-o. (41) 

Solutions of this equation must not diverge faster than the finite power of r and finite at 
r = O.The solution which satisfies this latter condition is derived by Whittaker equation 

(PW , 1 A 1/4 , 
for = l/2,whose solution is expressed by linear combination of M\^^|2 and Wa,i/2 

|x| \x\ 

S{x) is expanded near origin 

^WN^°= Mr(r:v2) +^'^'^- '''' 
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Therefore for normalizability of S{x) it is sufficent that l/r(l + A/2) vanishes. This condition 
is reahzed for 

1 + A/2 = = 0,1,2,..). (45) 

This is an eigenvalue condition for the couphng constant A. Thus fermion-antifermion(fa) 
system has discrete spectrum(A < 0).ln the case of positive copuhng the function M diverges 
as exp(x/2)[9].Thus fermion-antifermion(ff) system has no bound states(A > 0). Since the 
couphng constant has a dimension of mass,the mass of the ground sate is largest for fixed 
coupling A.Here we return to Ss{x) and SL{x).'We have 

CiM_„,i/2(2V2m|x|) , C2^-a,i/2(2V2m|x|) 
^S[X) = r-^ \ 7—: , i4bj 

D,M,,,/2{V2m\x\) D2Wa,,/2{V2m\x\) _ V2\' 

"^^^"^^ " R + R ' " - 

Eigenvalues are given as 

For the large distance in the case of these negative A"^, Ma,i/2 {V^rn \x\) blows up. Thus we 
choose Di = O.Here we have an approximate solution of the integral equation for short and 
long distance 

CiM_a,,/2{2V2m\x\) ^ 
S{x) ~ ^-^-j— (m \x\ <C 1) (48) 



\x\ 



D2Wa,y2{V2m\x\) 



\x\ 



{l<^m\x\). (49) 



We choose Ci = 1/2 for normalization of S{x).For small n we have an exphcit form of M, W 
in terms ol z — -\/2m \x\ 

Mi,i/2(^)/^ = exp(-^)2, 

M2,i/2{z)/z = exp(-z)2(l - z), {z < 1). (50) 

W^_i,i/2 {z)lz ~ exp(-z)/2z^ 
W^2M2{z)lz ~ exp(-^)/^3, (1 < z). 

Since at long distance the function W strongly dumps,we cut-off long distance and the 
solution may be approximated 

5(^)„ = MA,i/2(2|z|)/2|z|,A = n. (51) 
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FIG. 2: xf(x)(x),x^(x)(n) for m = 1 



If we transform into momentum space we get for n = 1,2, ^/2m = 1 



1 , 27 ■ p 



^^^^ ) = y ^ i^^(i^i)^ = wTw^' 'i-P 

, 167r(-l+p^ + 2p-7) 

) = ■ (52) 

Here we notice that the vector part of xi is not nor malizable, therefore we avoid n = 1 case 
for scalar-vector fa systems. Thus the ground state must be given hj n = 2. Correct solution 
in the whole region may be evaluated by numerical analysis of integral equation (32) with 
angular intergral 



m? — {2ttYp Jq \p ~ qI ' 4m 

with an input of Fourier transform of eq (49) 



in the right hand side of eq (53). We see the profile of X2{x) in FIG2. 



III. WARD-TAKAHASHI-IDENTITY FOR AXIALVECTOR CURRENTS 



In this section we examine the Ward-Takahashi-identity for the axialvector currents. In 
our approximation to BS equation fermion mass is assumed to be dynamical. Therefore if 
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we have a solution of the BS equation for psedoscalar it must satisfy the Ward-Takahashi- 
identity for the conservation of axialvector currents [13] 

limq^Spip')T,^ip',p)SFip) = 75} 7^ 0, (55) 

where Spip) is a solution of the Dyson- Schwinger equation or non-perturbative solution 
of the fermion propagator. The vertex function Fs^ has a massless pole g^/Q'^x''^(g), where 
X^(q') is a BS amplitude of psedoscalar Goldstone boson as a consequence of chiral sym- 
metry breaking. Following the notation in ref[5],f/(2) chiral symmetry is generated by 
{1,74,75,745} which is broken by dynamical fermion mass to a t/(l) x U{1) symmetry gen- 
erated by {/,745} for the degree of freedom of scalar and psedoscalar as {c, tt} .The set of 
7 matrices is 

' (73 \ I (Ti,2 \ I I 

7o = I h 71,2 




-aa 

745 = -«7475, {lix, 74} = 0, {7^, 75} = (56) 

Currents {7^4,7/^5} are a analog of axialvector currents and have a doublet of Goldstone 
boson which is called axial-scalar, 

Scalar-PS (745) BS amplitude is 

X'{P, q)^S + j-qV^ + j-PV' + e^^pP'^q'^Y^'T. (58) 

X^^(P, q) = S{P, g)745 + ?'^7,.45^' + P^l^^V + e.^pP'^q^YT. (59) 
To check the Ward-Takahashi-identity we consider the equation for axial-scalar 

im' + p')x^'{p) = ^J ciV^^^, (60) 



(27r)3y ^(p_p')2' 

AS { \ 

(m^ - U,)x^\x) = A^^^^-P, A^^ = A. (61) 

for ff system. Sign of the coupling is opposite to the one in even dimension.This is the same 
equation discussed in scalar-vector case and we have solutions 



Mx/2m,l/2' 


{2m \ 






\x\ 
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for A/2m = n(l, 2...). Following the Ward-Takahashi-identity x'^^i^) for is derived from 
the scalar part of the fermion propagator. Dy son- Schwinger equation for fermion propagator 
is written 

^^^^ = (2^ i ^ k^ + m- ' ^^^^ 

in the Landau gauge with linear approximation m — M(0).We define F{x) as the Fourier 
transformation of M{p) / [p^ + m^) 

^ / T^rXs exp^P • 2;)^— — ^, 65) 



2e^ / T^e^f'-^- 



d^P' iv'-x M{p') f d^k e/ 



(27r)3 p'^ + m^J (27r)3 F 

where k — p' — p and used that J d^p — J d^k.We arrive at a Schrodinger hke equation for 
F(x) 

(-a + m')F{x) = A^, A = ^. (67) 

Its solution is given as the groudstate of hydrogen atom 

2 3 
F{x) = ^exp(-m|x|),F(p) = j^^:^- (68) 

In this case the ground state solution corresponds to A/2m = l(m = e^/47r),where mass 
is largest for fixed coupling[10].We also meet this condition in the analysis of the fermion 
propagator based on low-energy theorem[ll,12].In fact we can determine a precise infrared 
behaviour of the propagator in the above analysis but do not have a definite ultraviolet 
behaviour. However if we demand the nonvanishment of the order parameter {ipip) 7^ 
,we obtain the condition for the anomalous dimension of the wave function which must 
be unity.This condition is approximately satisfied in the case of correction for photon 
propagator with vanishing bare mass in the Landau gauge[12].This is just the statement of 
Nambu-Goldstone theorem[14,15]. 
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IV. SUMMARY 



In this work we studied the solutions of spinor-spinor Bethe-Salpeter equation for massless 
boson in QED3 with ladder approximation in the Landau gauge.It is not easy to solve this 
equation directly in three dimension. However we use the approximate form of the integral 
equation for large and small momentum respectively. These are converted to the differential 
equation which are similar to the Schrodinger type equation. In this case an eigenvalue are 
determined by boundary condition at x — O.Thus we obtain solutions for massless scalar- 
vector fa systems with discrete spectrum.In our approximation massless scalar-vector ff 
system does not seem to exist in three dimension. Finally it is shown that Ward-Takahashi- 
identity for axial currents is satisfied with the Dyson-Schwinger equation for the fermion 
propagator in the linear approximation to dynamical mass . 
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